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A self-consistent-field theory is given for the electronic collective modes of a chain containing a 
finite number, N, of Coulomb-coupled spherical two-dimensional electron gases (S2DEs) arranged 
with their centers along a straight line, simulating a narrow micro-ribbon of metallic shells. The 
separation between nearest-neighbor shells is arbitrary and because of the quantization of the elec¬ 
tron energy levels due to their confinement to the spherical surface, all angular momenta L of the 
Coulomb excitations and their projections M on the quantization axis are coupled. However, for 
incoming light with a specific polarization, only one angular momentum quantum number is chosen. 
We show that when N = 3 the next-nearest-neighbor Coulomb coupling is larger than its value 
if they are located at opposite ends of a right-angle triangle forming the triad. Additionally, the 
frequencies of the plasma excitations depend on the orientation of the line joining them with respect 
to the axis of quantization since the magnetic field generated from the induced oscillating electric 
dipole moment on one sphere can couple to the induced magnetic dipole moment on another. 

PACS numbers: 73.20.-r, 73.20.Mf, 78.20.Bh, 78.67.Bf 


I. INTRODUCTION 

Over the years, plasmon excitations have been investigated both experimentally and theoretically for topologically 
different nanostructures such as the two-dimensional electron gas (2DEG) and for systems with layers of 2DEG 
forming a layered electron gas (LEG) systems [HE]. The LEG may be obtained in a multilayer semiconductor system 
by epitaxial growth such as GaAs/AlGaAs or in a type-II system such as GaSb AlSb InAs quantum-well structures 
containing an electron and a hole gas layer in each unit cell [3] . Another class of nanostructures for which the Goulomb 
excitations have been studied is the spherical 2DEG (S2DEG) |4]-[^ which has been employed in modeling metallic 
dimers laiH], clusters of carbon nanoparticles [9] and metallic chains of gold nanoparticles m^- However, since 
the spherical geometry allows for possible anisotropic coupling (unlike the cylindrical geometry [12] [13]) with an 
external light source by employing far-held polarization spectroscopy m, the Goulomb excitations for the S2DEG 
allow for more variation in frequency than the 2DEG by employing shells of different radii, forming hnite length chains 
with adjustable nearest-neighbor separations, or using bundles arranged in assorted conhgurations. Furthermore, we 
may model the S2DEG as an incompressible huid [8] or as a point dipole [14] in which only nearest-neighbor dipole 
interactions through an EM held are assumed dominant [14]. But, in such a model, the dipoles have been assumed 
to point in a chosen direction (like in an antenna array) so that only one component of the angular momentum is 
employed, e.g., if L = 1 then only one of M = 0, ±1 plays a role. In this paper, we overcome this restriction in 
investigating longitudinal and transverse plasmon-polariton modes in a chain with a hnite number N of metallic 
shells. Although our formalism may be applied to arbitrary N, we present results when N = 3 since it ahords us with 
the opportunity of comparing with the conhguration when the shells are positioned with their centers at the vertices 
of a right-angle triangle. This example demonstrates the orientational dependence of the Coulomb matrix elements 
as well as the ehects arising from the hnite length of a linear chain. 

Our model described below is drastically diherent from that in Ref. m because we include dynamically screened 
long-range Coulomb interactions and the orbital angular momenta and their projections onto the quantization axis 
are coupled. This means that the localized plasmons on each shell are coupled by the inter-sphere electron-electron 
Coulomb interaction. Also, an important diherence between the present model and the point-dipole model is that the 
electromagnetic held generated by time-varying charge distributions or current is retarded, whereas in our model the 
inter-shell coupling is electrostatic. Since the polarization function for the S2DEG vanishes for L = 0, there are no 
plasmon excitations for L = 0 [4]. This is in contrast with the model in Ref. m which allows plasmon excitations 
when the dipoles have a hxed direction, i.e., which corresponds to a chosen projection for an unspecihed angular 
momentum. 

With the use of our formalism, we may directly simulate the modes of oscillation observable by employing circularly- 
polarized light or a helical light beam. That is, the plasmons are detected experimentally using an external polarizing 
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FIG. 1: (Color online) Schematic illustration of a linear triad of S2DEGs. The angular momentum is L and its component is 
M. The shells are aligned along the x-axis. We consider two different cases of the quantization axis direction - along the x-axis 
(plot (a)) and along z-axis (plot (b)). The middle shell is centered in the origin, the distances between the centers are ai and 
a 2 , as indicated. 


field, which specifies a preferred direction, which is the axis in our case. Our primary goal is to consider a closely 
packed finite-length assembly of S2DEGs with all their centers on a straight line relative to the quantization axis. 
In such a system, the strength of the interaction between different nano shells varies with their separation as well as 
the relative orientation of the line joining their centers. With this in mind, it makes sense to distinguish the present 
investigation from both cases of a triad arranged at the vertices of a right-angle triangle [9] or an infinite linear array 
of periodically placed S2DEGs m- The assembly under consideration is shown schematically in Fig. We refer to 
each shell as j = 1, 2, 3, • • • , N, starting with the one on the far left. For N = 3, the middle shell is centered at the 
origin. We specify the distances between the S2DEGs as ai 2 = ai and a 23 = a 2 , so that ais = ai + a 2 . In all our 
numerical considerations, we investigate separately the cases corresponding to equal and unequal distances ai and a 2 
in conjunction with equal and unequal radii of the shells. We consider either the case of three identical shells with 
Ri = R 2 = Rs or the situation when the middle shell is larger than the other two: R 2 > Ri = R 3 . 


The spatial anisotropy of the plasmon excitations is just another example of the orientational dependence of the 
optical and transport properties of condensed matter systems. We now cite a few examples when such an anisotropy 
has been demonstrated. These include the electrical, thermal, mechanical and chemical properties of graphite along 
the a, b and c directions [15] and the elastic properties of carbon nanotube bundles [T6|. Also, the dispersion relation 
of the high-frequency optical 7r-plasmons for graphite was calculated by Chiu, et ah, m and it was shown that the 
plasmon excitations depend on whether the momentum transfer is parallel or perpendicular to the hexagonal plane 
lying in the Brillouin zone. The anisotropic conductivity of epitaxial graphene on SiC was reported in Ref. [18]. We 
note that there have been reports in which some of these anisotropic properties have been used in device applications. 
For example, in Ref. m, the tuning of surface plasmon frequencies to more efficient optical sensors was investigated. 


The rest of this paper is presented in the following way. In Sec. |TI] the theoretical formulation for deriving the 
plasmon equation for three spherical shells whose centers are on a straight line. The Coulomb matrix elements are 
shown to depend on the relative orientation of this line-of-centers with respect to the axis of quantization for the 
angular momentum associated with the spherical geometry. The plasmon equation in general involves the coupling 
of all angular momenta and their projections on the quantization axis. But, in Sec. |III[ we obtain the dielectric 


matrix assuming the impinging light sourse probing the Coulomb excitations has a specified polarization. Section IV 


is devoted to a presentation and discussion of our numerical results. We conclude with a summary of our results in 
Sec.jVj 


II. THEORETICAL FORMULATION OF THE PROBLEM 

We turn our attention to a system of N spherical shells with their centers on the x axis. The center of one of the 
spheres is located at x = 0 with radius Rq whereas the other spheres are centered at x =- ai and its radius is 
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Ri^ X = a2 and radius i? 2 , and so on. We assume no overlap of the shells so that we demand, for example, that the 
inequalities ai > i?i + Rq and a 2 > ao + a 2 are satisfied. In the absence of tunneling between the shells, the wave 
function for an electron on the j-th shell (j = 1, 2, 3, • • • , N) is given by 


< r I jV >= (f - (j - l)a 4 ), , 

fifi, 


( 1 ) 


with V = {l,m} and fj{r) = 5{r — Rj). The energy spectrum has the form 


_ h^l{l + 1 ) 

~ 2m* R] 


( 2 ) 


The equation of motion of the density matrix operator is 


= (3) 

where H = ~ is the Hamiltonian of the electron on the surface of the sphere, is the free electron Hamiltonian 

and ^ is the induced potential. The potential ^ satisfies Poisson’s equation 


V^^(r,cc;) = - Sn{T,uj) , (4) 

where and 55 is the uniform background dielectric constant. Additionally, (5n(r),co’ is the induced electron 

density. We use linear response theory to calculate the induced particle density as 


6 n{r,Lo) = '^'^ <r \ jiy >< jv \ pi(r,w) | j'v' >< j'v' | r > , 


( 5 ) 


where 


< I ^1(1*, cj) I j'v' >= 2 e 


HlU ^j'u' 


<jv I 4 >(r,u;) I j'v' > , 


(6) 


in terms of the Fermi-Dirac distribution function /o(e) and we express the induced potential as 4 >(r) = 

q' 

where Vb is a normalization volume. Then Eq. becomes 


2 e /o(ejiy) —/o(ej/iy/) 1 • ^ ^ / 1 ^ 

Sn{r,uj) = —2^2^ fc, ■ , . -^- <r\jv><jv |r> 


Vq T 


yv' 




( 7 ) 


or by taking the Fourier transform with respect to r 


( 5 n(q, uj) 


ho ^ + ^j'u' — ^jv 

J,3 

X jV> . 


ju > 


(8) 


The matrix elements < jv \ | j'v' > with wave functions < r | jV > given by Eq. Q may be calculated using 

the expansion of a plane wave in spherical waves 
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e*'!- =47r^ 

l,m 


where ji{x) is a spherical Bessel function. The result is 


(9) 


< I I jV' >= 47r%/e*(^- 
X ^ i^jLiqRj)Ylj^icD [ an YCMYLMm Yi,^,{^l) , 

L,M 


Substituting Eq. (10) into Eq. we obtain after some algebra 


( 10 ) 


/o(ej,0 - fo{ej,i') 


Sn{q,uj) = {4Trf^ y^y^Yl ^ 

{-i)^jL{qRi)YLM{ii) [dn Yir^,{Q)YlM{Q)Yi^{n) 

T T\/r ^ 


L,M 






Y i^'3L'{q'Rj)Yl,M'{i() [ dn Yil,{n)YL,M'm Yi,^,{Q) . 

L',M' ^ 


( 11 ) 


(5n(q, cj) = (47r)^^ 

^%tir 2 rv 


/o(ei,0 /o(ei,r) (2; + i)(2r + l) ('J ^ 


X jL(gi?,)UM(q) 

M 

X E (4, g;, q,) ji(g'i?,)UV(q ) • 


( 12 ) 


in terms of Clebsch-Gordon coefficients. Taking the Eourier transform of Eq. (|^ we have 4>(q) = —4:7reSn{q)/s. 
Using this relation in Eq. (11), we obtain 




(5n(q, cj) = — 


87re^ 


(47r)2 ^ n,-i(w)e-'(l-i)«*“ jL(5i?,)[/j,L,MnM(q) 

j,L,M 


(13) 


where Uj^Li^) is the density response function of the j-th nano shell with 
where 


liLfCJj = > - --- 

^^ huj T 6// — 6/ 
1 , 1 ' ^ ^ 


(14) 


and 


a 


j,LM 


y.x ' ”1” 


LxLyLz 


(15) 


qx,qy,qz 


Here, L^^Ly and are normalization lengths, with Vb = L^LyL^. Substituting the expression for Sn{q) given in 


Eq. (13) into Eq. (15), we obtain 
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FIG. 2: (Color online) Coulomb matrix elements V 12 (in units of 2e^/{'KegR) — 5.7 el/ ) for a pair of interacting S2DEGs as a 
function of the separation a (in units of their radius = 1 nm) between their centers. The displayed matrix elements are Vbo,oo 
(red curves), Vio,io (green curves), (blue curves), and (black curves). Plot (a) shows the case of x—alignment, and 

plot (5) shows results for z—alignment. The largest matrix element, Vbo,oo, which corresponds to L = L' = 0 and M = M' = 0 
does not contribute to the plasmon excitation spectrum. However, it is given for comparison. The insets demonstrate the 
corresponding interaction potential matrix elements (same color) as a function of the ratio R 2 /Ri of the radii of the two shells 
for fixed separation a between their centers. 


N 00 L' 

EE E 

j' = l L'=0 M' = -L' 


2e^ 


Uj'^L'M' = 0 , 


(16) 


where 


We have explicitly by setting j = 1, 2 or j = 3 in turn for each of the three spheres 


(17) 


1 


2e 


2e^ 


£g{2L T 1)77- 


n 


j,LM 


(uj) 


U: 


j,LM 


- T/ E ^j',L'{f^)Vj'L'M',jLM{Rj,Rj',0')Ujr^L'M' — 0, 




(18) 


If the spheres are identical, then we need only set j = 1, but still have to do the sum over j' = 1,2,3. The set 
of linear equations (18) have nontrivial solutions provided the determinant of the coefficient matrix of {Uj^lm} is 
zero. Consequently, ^smon modes with different values of L on the two shells may now be coupled via the Coulomb 
interaction. Since Vilm,L'M'{Ri^ R 2 ', cl) ^ 0 in the limit a ^ oc, this matrix is diagonal when a ^ i7i,i72 and the 
plasmon mode equation reduces to the result for isolated shells 




( 19 ) 
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FIG. 3: (Color online) Particle-hole modes and plasmon excitations for a single S2DEG and for a linear assembly of three 
shells as a function of their radii. Panel (a) shows the particle-hole modes for a single S2DEG. Panel (b) presents the plasmon 
excitation frequency of a single S2DEG as a function of its radius. All the other plots in panels (c) through (/) demonstrate the 
dependence of the plasmon frequency as a function of radius R for three equivalent S2DEGs with Ri = R 2 = R and with equal 
separations ai = a 2 = a = 2R + 0.1 nm. Panels (c) and (d) represent the plasmon branches of an x—aligned linear triad with 
various ranges of their radii, i.e., from 1.0 to 1.7 nm for (c) and 1.0 to 2.4 nm for (d). Plots in (e) and (/) show the corresponding 
plasmon dependence for the case of 2 :—alignment. While plot (e) shows results for radii in the range from 1.1 to 1.5 nm, and 
plot (/) gives the dependence over a narrower range 1.2 to 1.3 nm. 


where is the dielectric function for the j-th shell. The significance of equations (18) for chosen L,M is that 

they give explicitly the effect of the Coulomb interaction on each shell through as well as the coupling between the 
pair of shells through the Coulomb matrix elements Vj'L'M'jLMiRj, Rj' ^ct). Additionally, the nature of this coupling 
may be characterized in the following way when carrying out numerical calculations. For chosen L and M satisfying 
—L < M < L, there are 2(2L + 1) x 2(21/ + 1 ) x 2(2L + 1) elements in a block sub-matrix which includes (2L + 1) 
elements along the diagonal, equal to (21/ + 1) diagonal elements equal to £ 2 ,l{^) and {2L + 1) diagonal 

elements equal to For example, if we consider the coupling between sub-matrices with angular momentum 

L = 1, 2, 3, • • • , A/", then the dimension of the matrix is 3 J])^^^(2I/ + 1 ) = 3N^ + 6N. Specifically, if we use just the 
L = 1 sub-matrix, we have a 9 x 9 matrix which we discuss below. 
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III. DIELECTRIC FUNCTION MATRIX CALCULATION 


We now restrict our attention to the case when N = 3 and the angular momentum quantum number is chosen as 
L = 1. The zeros of the dielectric function for the three aligned S2DEGs, which determine the plasmon excitation 

•f-)- 

energies of the system, may be obtained as the solution of a determinantal equation Det TWg = 0, where is a 
complete 9x9 matrix, including the dielectric function of each nanosphere, as well as the inter-sphere interaction 
potential elements. In our formalism above, we took the spherical shells as being centered on the x-axis. However, 
without loss of generality, the spheres may be placed along the axis of quantization, the z-axis. In either case, the 
dielectric matrix may be presented as a block matrix as follows: 




M, 




Vf3\ 



^2 


(20) 

vfi 


/ 



where, in our notation, the submatrix V- is diagonal. Each diagonal block consists of three identical elements, 
representing the dielectric function of an isolated S2DEG: We have 




^L=i{Ri, cj) 0 0 

0 eL=i{Ri, Uj) 0 

0 0 eL=i{Ri, uj) 


( 21 ) 


for i = 1, 2 or 3. The coupling between different shells is given by the off-diagonal blocks Vj^-, and representing 
the interaction potentials between the adjacent (1 — 2) and (2 — 3) shells, and next-nearest-neighbors (1 — 3). This 
geometrical arrangement leads to the main specific feature of our case of the aligned shells. In this regard, we consider 
an off-diagonal 3x3 interaction matrix Yjj ', which may denote either Yfj or Yjj and include their common features 
as demonstrated in Eig. We emphasize that the analytical formula of the matrix and the expression for each 
elements are identical for a next-nearest-neighbor or an adjacent pair of S2DEG’s, while the only difference comes 
from the distance between the shells, resulting in the value ai 2 + a 23 for the next-nearest-neighbor pair. Specifically, 
for L = L' = 1, the matrix consists of nine elements determined by the M and M' as follows: 


V: 


MM' 

33' 


{R]ajj>) = IVl=i{R, 


kM=-l,M' = -l kM=-l,M'=0 = l 

kM=0,M' = -l Vm=0,M'=0 Vm=0,M' = 1 

Vm=i,m'=-i Vm=i,m'=o Vm=i,m'=i 


( 22 ) 


Each matrix elements depends on the orientation of the nanoshells and will be calculated below according to Eq. 0- 


Each matrix element in Eq. (21) may be expressed as eL=i(i?i, cj) = 1 + 2e^/{3esR)IlL=i{^)- Consequently, for 


the simplest case of three non-interacting shells of equal radius R, the plasmon branches are obtained by solving for 
the zeros of the function 


Bet M, = {Bet Vf = (1 + 2e‘^/{3e,R) nL=i(w))^ (23) 

i.e. we obtain a ninefold degenerate plasmon solution for a single isolated S2DEG. This solution corresponds to the 
result for interacting spheres as a function of their separation when both a\ 2 lRi cxo and a 23 lRi ^ oc. In proceeding 
to take the inter-shell Coulomb coupling into account, we must distinguish between the and x—alignments, leading 
to different analytical forms for the interaction submatrix. Let us now consider the interaction matrix elements for both 
cases of alignment. We conclude that the result of the integration is non-zero only for certain elements, representing 
an interesting set of selection rules. The situations varies for each alignment. Eor x—alignment, the order of the 
remaining Bessel function is determined by the difference \M — M'\ and the elements are classified correspondingly 
(see Eq. (HI))- Eor example, all the elements with \M — M'\ = 1 are equal to zero due to the inherent symmetry of 
the polar integral. Eor 2 :—alignment, the azimuthal angle (j)— is determined by the phases of the spherical harmonics 
in Eq. ( p!7[ ), so that only elements with M = M' are non-zero, i.e. we obtain only two such elements. 

We first consider the simplest case of three identical shells each of radius R and with the same number of filled 
energy levels or Eermi level Lp- We set the separation ai 2 = a 23 = a, so that ai 3 = 2a. In previous work, we obtained 
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FIG. 4: (Color online) Plasmon excitations for a linear assembly of three identical S2DEGs as a function of the separation a 
(ai = a 2 = a). The radius of the spheres was chosen to be i? = i?i = i ?2 = 0.8nm. The upper panels (a) and (b) correspond 
to the case of x—alignment, and the lower ones (c) and (d) to z—alignment. Each panel presents results for a specific range of 
the separation a: 1.7to5.0nm for (a) and (5), respectively. In panel (c), the range is 1.6tol.9nm and l.GtoG.Onm in panel 
(d). The solution for a system of non-interacting shells Qo = 9.17 eV/h is indicated for comparison. 


[7] the pair-wise interaction potential matrix for two S2GEGS with radii Ri and 7^2, aligned along the z-axis with 
separation a 


f VI 0 0\ 

Vjjv = nL=i(i?, w) 0 i;3 0 , 

\ 0 0 Vi J 

where 


(24) 


vi = Vipii 


'^3 — ^ 10,10 


6e^ 

ires 

6e^ 

Tre^ 



{sm{aq) - aqcos{aq)} ji{Riq)ji{R2q) 


OC 

/ djQ 

—p- \^2aqcos{aq) + ((<a^)^ — l) sm{a q)} 

0 


(25) 


The Of—dependence accounts for the fact that two adjacent spheres (1 — 2 and 2 — 3) have equal separation a, whereas 
the far-removed pair 1 and 3 have their centers separated by distance 2a. 

For a ^—aligned triad, each interaction matrix block matrix consists of only two different Coulomb potential matrix 
elements and, most crucially, is diagonal. This established fact results in the factorization of the determinant of the 
A4e matrix. The resulting product matrix is 


4 

Det Me=Wmx 

A=1 


(26) 


with 
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FIG. 5: (Color online) Exact numerical solutions for the plasmon frequency for three z—aligned identical S2DEGs with equal 
radii Ri = R 2 = R and separation between them m = a 2 = a. The plasmon frequency dependence on the radius of the spheres 
for the range 1.3nm < R < lAnm is shown at plot (a). The separation between their centers is a = 2R0.1 nm. Plot 
(b) represents the frequencies of the plasmon branches for the three spheres with R = 0.8 nm as a function of the separation 
a. The solution for a system of non-interacting shells = 9.17 eV/h is also provided for reference. Since the determinantal 
equation may be factorized, we clearly separate for different solutions 1 — 4, although two of them are doubly degenerate (1 
and 3), and each of the other two solutions 3 and 4 yields two plasmon branches (four in total). 


/ 2e^ 

mi = I eL=i{R, oj) - UL=i{u})vi{R,2a) 

V 

2e2 

m 2 = eL=i{R, w)- IlL=i{u)vs{R,2a) 

( 2e^ /2e^ 

m 3 = uj) + eL=i(i?, to) -nL=i(w)i;i(i?, 2 a) -2 ( - IlL=iVi{a,R) 


2e 


m 4 = ei=i{R, oj) + eL=i{R, to) -nL=i(a;)v 3 (i?, 2 a) - 2 - IiL=i'V3{a, R) 


2e^ 


(27) 

(28) 

(29) 

(30) 


each “factored linear” multiplier results in one solution, either singly or doubly degenerate, whereas each “quadratic” 
part gives rise to two different solutions. 

For X—aligned shells, the potential submatrix is not diagonal, so similar factorization is not possible and we have 


Vi 0 V2 

^jj'=^ l=i{R, oj) \ 0 i;3 0 

V2 0 Vi 


(31) 


with the elements as follows: 


vi = Vi_i,i_i = 


3e^ f djQ^ 

Vii,ii = — J {aqcos{aq) + {{aqf - l) sm{aq)} ji{Riq)ji{R 2 q) 

0 

C» 

3g^ r (f(^ 

V 2 = = - / -—^ {3aqcos{aq) + {{aqf - 3) sin{aq)] ji{Riq)ji{R 2 q) 

0 

oo 

f djQ 

V3 = Uo.io =- / 7 — 7 ^ {aqcos{aq) + sin(ag)} ji(i?ig)ji(i? 2 g) 

nCa J [aqR 


(32) 


(33) 


As the final step of our discussion of the potential elements, we address the properties of the largest one, corre¬ 
sponding to 1/ = L' = 0 and M = M' = 0. It does not contribute to the plasmons of our system because ni;,=o(^^) = 0- 
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Clearly, it demonstrates complete spherical symmetry, i.e. this is the only element, having this same property for both 
X— and z— alignments. Also, due to its specific symmetry properties, this is the largest element. Some straightforward 
integration in Eq. 18 leads to the following expression: 


Vbo,oo(^i, ^2, ci) 


3e^ 

TTCs 


oc 

j dqjo{Riq)jo{R2q)jo{aq) 

0 


(34) 


In the following section, we employ these results to determine numerically the plasmon excitations when three shells 
are lined up along or perpendicular to the axis of quantization. 


IV. NUMERICAL RESULTS 


As the first step for our numerical results, we investigate the Coulomb interaction potential matrix elements 
Vj'^L',M 'which are involved in determining the plasmon excitations of the assembly according 
to Eq. ( [l7| ). The calculation of the potential matrix elements of each interacting pair of S2DEGs is similar to our 
previous considerations HIE]. All non-zero potential matrix elements are given in Eig.j^for L = 1. We note the each 
inter-sphere matrix element has a strong dependence on the distance between the shells, so that the interaction is 
significant for a closely packed assembly with ajj' ^ Rj + Rj >. This behavior is similar to the single shell potential 
matrix element ^ 1/[{2L + ^)Rj]- The magnitude of each element at a given distance depends on the symmetry of 
the quantum state, i.e. the values of L, L', as well as M and M'. The largest one corresponds to rotational symmetry 
when L = L' = 0 and M = M' = 0. However, it does not couple with irradiation with specific polarization and, there¬ 
fore does not contribute to the plasmon excitation since the polarization function is identically zero for zero angular 
momentum. All other elements are of the same order of magnitude. Existence of non-zero potential matrix elements 
are given by specific selection rules, representing an interesting contribution of our work. Thus, due to symmetry of 
the interaction, the only non-zero elements for z-alignment are those with M = M', which is not the case for the 
alternative situation of x—alignment. Only the elements with L = L' = 1 contribute to the plasmon spectrum for both 
cases of shell alignment. As a result, we are left with three different non-zero interaction matrix elements in the case 
of X—alignment, while the shells aligned along the 2 : axis exhibit only two of them, which is reflected in constructing 
the interaction matrices in Eq. (20). The insets in Eig. [^show the behavior of each potential matrix element for the 
case of unequal radii of the interacting shells. We conclude that there is no pattern for such dependence, except that 
each element, either positive or negative, increases (the absolute value) as a function of p = ^ 2/^1 for approximately 
equal radii Ri ^ ^ 2 - 


We now turn to a discussion of our numerical results for the plasmons of a linear triad of S 2 DEGs for angular 
momentum L = 1. The plasmon excitation energies of three Goulomb-coupled S 2 DEGs are obtained by solving for 
the zeros of the determinant of the matrix in Eq. (26) when the S 2 DEGs are centered on the z axis and solving the 


general equation (20) with the interaction submatrix given in Eq. 


[M| ). In general, we must deal with a matrix of 
18l since all angular momenta corresponding to 


infinite order whose elements are given by the coefficients of Eq. 

L = 1, 2 • • • , 00 are coupled. However, for incoming light with a specific polarization, only one angular momentum 
quantum number might be needed in carrying out the calculation to determine the plasmon frequencies since in 
general the angular momentum of light may be carried by either orbital motion (helicity) or spin motion (circular 
polarization). If a finite angular momentum of light with L = 1 is used for incidence, the magnetic field generated 
from the induced oscillating electric dipole moment on one sphere can couple to the induced magnetic dipole moment 
on another displaced sphere. This results in a 9 x 9 determinantal equation with elements corresponding to j = 1, 2 
and 3 in Eq. (18) with L = 1 and M = 0 ±1. The resulting matrix consists of thee 3x3 diagonal blocks, related 
to the plasmons of each individual S 2 DEG and the remaining six off-diagonal 3x3 submatrices with j 7 ^ j' = 1 , 2 
and 3. The diagonal block could be either identical for the case of three equivalent shells, or different, depending on 
the radius and the number of electrons on each S2DEG. Each block is diagonal, which means that for non-interacting 
shells the total matrix is diagonal and its determinant may be easily factorized. 


Both the intra- and inter-shell Goulomb interactions contribute to the plasmon excitation energies. Eirst, we 
determine how the plasma frequencies depend on the radius of three identical S 2 DEGs. Our results are presented in 
Eig. i Eirst, we present the single-particle excitation region (SPER) in the uu-R plane as being identified by finite 
imaginary part of the non-interacting electron polarization function with the R—dependence arising from 

the single-particle energies in Eq. §. The SPER, which is also referred to as the particle-hole excitation region, 
shows the u-R regions where there is natural damping of the plasmon modes. We see by comparing our plots in Eigs. 
[^a) andj^b) that the plasmons for a single S2DEG are undamped, well separated from the SPER and their energies 
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decrease monotonically with increasing radius. We present our results in Figs.j^c) through [^f) for interacting shells. 
The two distinct cases of either x— or 2;—alignment clearly lead to different solutions for the plasmon modes. Whereas 
for X—alignment, there are nine branches, there are only six solutions for the plasmons in a ^—aligned assembly. The 
corresponding a—dependence (vs. the distance between the S2DEG centers) is displayed in Figs.j^a) through [^d). 
In general, the plasmon modes do not have uniform intensity. Specifically, the peaks are not strong for the lower 
frequency modes because of their proximity to the particle-hole mode region. We found that when the separation a 
increases, the plasmon branches tend to the single degenerate solution denoted by flo? corresponding to non-interacting 
shells. 


The mathematically simpler case when there is 2:—alignment, which allows one to factorize the determinantal 
equation for the matrix in Eq. (26), deserves some attention. Our calculations yield each plasmon mode as a numerical 
solution of the determinantal equation, and present our results in Eig. Once again, we confirm that there are six 
distinct branches, corresponding to the four factored matrices, two of them are quadratic resulting in two solutions 
each. Two of the solutions are doubly-degenerate. The branches are not symmetric in intensity with respect to the 
horizontal line flo- 

So far, in this paper, we presented numerical results when the radii of the spherical shells are equal and also when 
the separations between nearest-neighbors are equal. Figures [^a) through [^d) examine the behavior of the plasmon 
modes of a linear assembly with unequal nearest-neighbor separations as well as when the S2DEGs have different 
radii and chemical potentials, which corresponds to more realistic situation of fullerene aggregates. First, in this 
regard, we have investigated the mathematically simpler case of z—aligned shells with equal radius i? = 0.8 nm but 
unequal distances between their centers. In this regard, there is no symmetry of the interaction submatrices and 
we may no longer use the matrix factorization in Eq. (26). Instead, we obtain a lower- order factorization, which 


consists of two multipliers, each of cubic order. Therefore, we still have six solutions, but there is no degeneracy of 
the frequencies. Once the ratio 7 = a2/ai is increased, the electrostatic interaction between the nearest-neighbor 
shells which we refer to as (2) and (3), and naturally between next-nearest-neighbor shells (I) and (3) is decreased, 
with the plasmon frequency of the assembly tending to the single-shell plasmon mode frequency Oq = 9.17, eE//i. 
The remaining plasmon branches are similar to those of the interaction between (I) and (2). These branches are 
symmetric with respect to the asymptotic line flo, as we obtained for the case of two interacting 2;-aligned shells [7]. 
When the distances between nearest-neighbor S2DEGs are comparable, the plasmon modes do not behave as 
which is a new feature compared to all the previously considered cases HIE]. 

Finally we address the case when the radius of the middle shell (2) is larger than the radius of the other two: 
R2 ^ R\ — 7^3 with p — ^ I* Now we have two asymptotic lines Hq ^ and Hq \ The larger solution Hq 

corresponds to the two spheres with lower radius Ri = i?3, i.e. is double degenerate, which results in a stronger peak 
(see Fig.[^(c) and (d)). The plasmon modes are strongly asymmetric. We also conclude that only two branches have 
their asymptotes along the line Hq , corresponding to the middle shell (2). 


V. CONCLUDING REMARKS 

In this paper, we carried out a model calculation of the plasmon relation for a narrow ribbon of fullerene atoms 
or metallic shells which were simulated by a linear array of S2DEGs. In neglecting the width of the ribbon, we only 
investigated effects arising from the edges along the length and not the width of the narrow ribbon. The coupling 
between plasmon excitations for the 2D electron gases on the surfaces of an finite number of shells is introduced using 
the random-phase approximation. We demonstrated that the quantization of the energy levels for each 3DEG on 
the surface of a spherical shell leads to Coulomb matrix elements depending on the angular momentum L and its 
component M. These Coulomb matrix elements represent the inter-shell coupling of the plasmon excitations and are 
clearly anisotropic, depending on the angle the line-of-centers makes with the quantization axis leading to anisotropy 
in the plasmon coupling with respect to the direction of the probe field. Furthermore, we proved that the strength 
of the inter-shell Coulomb interaction does not only depend on the distance between shells but also on the direction 
of the line joining these centers. We presented detailed results for the plasmon excitation frequencies as functions of 
the separation between shells and their radii. The effects of unequal separation between nearest-neighbor shells and 
when the shells have unequal radii were also investigated. 

From an experimental point of view, when light with a specified finite orbital or spin angular momentum is incident 
on the ribbon, the magnetic field generated from an induced oscillating electric dipole on any sphere may couple to 
an induced magnetic dipole on one of the spheres in the array where the coupling strength is determined by their 
orientational direction either parallel or perpendicular to the probe E field. This leads to dimerization of pairs of 
S2DEGs confined on two separate shells. Therefore, the spectra of the plasma excitations are different in the cases 
when the quantization axis is parallel or perpendicular to the array axis. 
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FIG. 6: (Color online). Plasmons for an asymmetric linear assembly of S2DEG’s. The upper panel demonstrates the plasmon 
solutions for three equivalent fullerenes with different distances between their centers: a 2 = 1.2 ai for plot (a) and a 2 = 1.2ai - 
for (6). Distances are specified according to the schematics in Fig[^ The range of the smaller distance ai = ai 2 is 1.6 — 2.6 nm 
for both plots (a) and (b). The solution for a system of non-interacting shells Qo = 9.17 eVjh, being the asymptotic solution, is 
provided for comparison. The lower panels (c) and (d) demonstrate the plasmons for three different spheres with R2 ^ Ri = i? 3 , 
7 ^ 2 / 7 ^! = 1-05 and 1.02, correspondingly. The distances between the spheres are equal: ai = a 2 = a. Each plot has two 
asymptotic solutions - = 9.17 eV/h, corresponding to Ri = 0.8 nm, = 8.48 eV/h for R 2 = 0.84 nm (panel (c)) and 

= 8.87 eV/h for R 2 = 0.82nm (panel (d)). The distance a between the centers of the spheres is displayed in the range of 
1.75 — 2.2b nm for both plots (a) and (b). We chose the case of z-alignment for all the plots. 


This work covers fundamental aspects such as anisotropy, many-particle quantum effects and electron-plasmon 
interactions for a novel plasmonic material. Applications such as biosensors for health care, devices for telecommu¬ 
nications, and near-field instrumentation may be explored. We note that our model calculations may not only be 
applicable to metallic particulates but also have broader implications to fullerenes. The numerical results we derived 
demonstrate significant new information in the area of plasmonics and are experimentally observable.. In this con¬ 
nection, there have already been some experimental measurements showing similar effects in nanoparticles [20]. Also, 
our work has a bearing on that of hybridization for surface plasmons in metallic dimers Hi in which the plasmon 
frequencies depend on whether the quantization axis is parallel or perpendicular to the inter-particle axis. 
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